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Abstract 

We introduce new methods for investigating propagation of electrons in a multi band system 
with spin orbit interaction which are time reversal invariant . Therefore Kramer's theorem imposes 
constraints which give rise to non-trivial Berry connections. As a result chiral zero modes appear 
at the interface between two domains characterized by parameters which are above or below some 
critical values. The variation of the parameters is due to disorder, geometry or topological disorder 
such as dislocations and disclinations. The mechanism might explain the high conductivity coming 
from the bulk of the topological insulators. We introduced the method of curved geometry and 
study the effect of dislocations on the Topological Insulators. As a demonstration of the emergent 
Majorana Fermions we consider a P-wave wire coupled to two metallic rings. 
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I-INTRODUCTION: THE BERRY CONNECTION AND CURVATURE 
FOR PARTICLES IN SOLIDS 



One of the important concepts in Condensed Matter Physics is the idea of Topological 



order [l4lOl. 



12 



32| . This concept has been confirmed experimentally by the 3D class of 
Bi2Te^ and Bii^^Sbx (characterized by a single Dirac 



Topological Insulators (TI) BioSe 



cone which lies in a gap) 



15 



18 



19|. At the boundary of the 3D Topological insulator 



one obtains a 2D surface which an odd number of chiral excitations. The two-dimensional 

CdTe/ HgTe/CdTe quantum wells behave as a 2D TI when the width of the well exceeds 

I 

6.3 nm, the s type Tq band to be below the p type Fg 30?]. The boundary of the 2D surface 
is described by a ID gap less chiral mode. (It is the single chirality which prohibits back- 
scattering and gives rise to the quantized spin-Hall effect.) Recently the QSH (Quantum 
Spin Hall effect) 30] in CdTe/ HgTe/CdTe has been discovered jisl, 31] and at the edges a 
spin Hall state dubbed helical liquid has been proposed. This spin model is time reversal 
invariant and therefore obeys Kramer's theorem = — 1 (T is the time reversal symmetry 
T = —ia'^K, K is the conjugation operation and cr^ is a Pauli spin operator). The helical 
liquid is given by a one dimensional massless chiral fermion: right fermions with spin up 
and left fermions with spin down. The spin Hall Hamiltonian is a physical realization of 
the lattice model proposed by 

In a variety of Physical problems these properties are described by the topological Chern 
number, Chn = jrpD=2nTr[-^{^Y] where T^=2" represents the D dimensional Torus ( the 
First Brillouin Zone (FEZ) )in D dimensions and F represents the curvature given in terms 
of the ground state eigen-spinors . 

In a physical language this means that ground state is protected by a gap which is 
stable to weak perturbations. At the boundary of such a system one finds an odd number 
of chiral massless fermions. This is revealed by the presence of a single Dirac cone. In the 
absence of external magnetic fields a well known theorem the no go theorem prohibits the 
existence of a single Dirac cone. The resolution of this paradox is due to Wilson js] which 
showed that by adding terms of the form '^f=i[cos{Kia) — 1] to the Dirac Hamiltonian 



one can choose the parameter of the model such that one has an odd number o 
fermions on the boundary. This situation has been realized in the class of materials 



chiral 



15 



Bi2Ses, Bi2Te^ and Bii^^Sbx- Theoretically a number of methods have been proposed for 
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identifying the Chern numbers js-S] and the Chern — Simons coefficients which describe 
the electromagnetic response of the system. The relation between the space dimension D 
and the index n is D = 2n. The ffist Chern number n = 1 exists for D = 2 and is non zero 
for system with broken time reversal symmetry (Quantum Hall effect) H]- The second 
Chern number Chn=2 is non zero for D = A and is not restricted to the case where time 
reversal symmetry is broken. The time reversal invariants system (TRI) in the presence 
of the spin orbit interaction can give rise to topological insulators TI. Recently \}\ has 
proposed a method for computing the topological invariants for non interacting electrons in 
periodic crystals which are time reversal invariant {TRI) and obey the inversion symmetry . 
For system with inversion symmetry the TRI are given in terms of the parity of each pair of 
the Kramers degenerate occupied bands (at the four (D=2) and eight (D=3) time -reversal 
and parity invariant points in the Brillouin zone). A general theory based on particle-hole 
symmetrj^time-r ever sal symmetry and their product allows to identify 10 symmetry class 
14l |. However we believe that the calculation of the invariants for systems which do not 
have parity invariance can be done, directly by computing the zero modes on the boundary. 
Using the index theorem we can relate the number of zero modes with the Chern numbers. 

Recent experiments point to interesting Quantum oscillations in topological insulators 33 1 
evidence of Landau Quantization in TI 3^ and Aharonov-Bohm oscillations fo r cy linders 



35|. Some of the experiments demand the use of geometry in curved space 

[36|,|3iQ- It 



seems that topological effects such as dislocations and disclinations might play an important 



role 
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Zero modes can appear also at the interface between Superconductors and normal metals. 
For example the Chiral superconductor {p^ + ipy) in the weak pairing limit is given by a 
similar hamiltonian as considered for the spin Hall effect. The physical difference is that 
for the superconductor the spinor is built from particles and holes and the negative 
energy of particles describes the same excitations as positive ones. As a result, 
the creation and annihilation operators are given by the Majorana fermions j4o| built from 
a linear superposition of a fermion with the charge conjugated fermion. The edge states 
for this case will be given by the one dimensional Ising model which has zero modes with a 
degenerate Z2 ground state 41j (Majorana zero mode is trapped in each vortex core which 
leads to a ground state degeneracy of 2""^ in the presence of 2n vortices). 
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II-THE METHOD FOR TOPOLOGICAL INSULATORS 



In this section we will introduce a set of methods which we will ue in this presentation. 
We will demonstrate our methodolo gy using a simplified version of the spin orbit model 



studied by 



Qflflfl 
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ITj . We consider a four dimensional lattice in which 



we introduce a defect in one of the direction. We show that for certain parameters a zero 
mode appears at the interface between domains which have different parameters. When 
we take into consideration the time reversal points in the Brillouin zone we can introduce 
an alternative way for determinig the criteria of topological insulator (TJ). In section B2 
we characterize the TI in terms of the Chern numbers. We show that the different Chern 
numbers are determined by the symmetry group which connect different different region in 
the Brillouin zone. In section B3 we show that the Hubbard Stratonovici field which results 
from the electron-electron interaction can be treated as fourth dimensions giving rise to a 
TI in three dimensions. 



Using the momentum space representation in periodic solids 23 



45 



471] we have con- 
structed a gauge theory in the momentum space which gives rise to non-commuting 



221] coordinates . 

For periodic potentials the eigenvalues are characterized by the Bloch eigenfunctions 



U 



n,K 



(q) and eigenvalues En{K) where n is the band index. In the momentum repre- 
sentation the coordinate r* is represented by idxi- When more than one atom per unit 
cell are involved we can have band crossing with linear dispersion (in the vicinity of the 
crossing points). Due to the band crossing, we can not use a single matrix connection 



< K\d\K >=< K\dKi\K > dK' { \K > are the eigenvectors) for the entire first Bril- 



louin zone 
connection 



i^BZ) (which is equivalent to the D dimensional torus T ). Due to the multiple 



9|,l20|, 



22 



24 



I, |28| one has to introduce transition functions between the 
different regions in the FBZ. When the number of bands crossing is odd we have an odd 
number of Dirac cones. For time reversal invariant (TRI) systems (with the the time re- 
versal (TR) operator T) we can apply Kramer's theorem = — 1 at the special time 
reversal points K = 0, K* = G — K* (where G is the reciprocal lattice vector). We find 
a number of materials which obeys Kramer's theorem with an odd number of Dirac points 
Bi2X3,{X = Se,Te) , Bii^^Shx coined Topological Insulators TI. These materials are char- 
acterized by strong spin-orbit scattering and have orthorombic crystal structure with the 



5 



space group D^^ and have a layer structure with five atomic layers per unit cell. Using a 



projected atomic Wannier function we can write a projected Hami 
sixporbitals \px,t>, \Py,t> , \PzA> and \Px,l>, \Py,i> \Pz,l> 



tonian in the basis of the 



3- 



Using these orbitals 



we obtain a six band Bloch Hamiltonian which are further projected to only four bands. 

IIA-The effective model and the emergent chiral modes at the interface be- 
tween two domains 

In order to demonstrate the physics involved we will consider the toy model of two orbital 
per unit cell : 



H = iXso^[ia]aj+i + a''j^-^aj) - {bpj+i + b]^^bj)]+M^{a]bj + blaj)-t^[{albj+i + a]^^bj) 

j 3 3 

(1) 

The eigenvalues are given by: e(i?) = ±^{2\soSinKY + (M - 2tcosKY. 
This model has the interesting property that for M — 2t = it develops a Dirac cone at 
= ( at M + 2t = the Dirac cone is at K = ±7r). For M - 2t < have at = a gap 
given by 2|M - 2t\. 

Next we extend the one dimensional model to model which has a boundary at x = 0. We 
assume that the Hamiltonian parameters obey for x < , (M — 2t)\x^Q < and for x > 
the parameter are (M — 2t)\^yQ > 0. As a result of the fact that the Hamiltonian for a; < 
is different from the Hamiltonian for x > gives rise to a mid gap state with zero energy 
which is to a zero mode. 

In order to demonstrate that numbers of zero mode is odd and are chiralic we will we 
will include a spin orbit interaction . We introduce the basis | t> ®|1 >, I i> ®|1 >) 
I t> ®|2 >, I t> ®|2 >• We will use the Pauli matrix r to describe the band and a to 
describe the spin. In this basis we can describe the spin orbit interaction by the term: 

D 

iXso ^■^(r)(a* ® ri)^(f + a^) + h.c. (2) 

i=l 

The energy difference between the two orbitals is represented by 

M^^(f)(/(g)ri)^(f) (3) 
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The hopping energy between the different orbitals is given by 



-t ^ ^^(r)(/ (g) ri)*(f + Oj) + /i.e. (4) 

i=l,2,..D 

This Hamiltonian captures the same physics as used by [71, [l2|, [l3| . The Hamiltonian can 
be rewritten using the Dirac matrices with 7" given in the chiral representation . We have: 
^1 = 7°(o"-'^ (g) r3),7^ = 7°(cr^ (g) Ta), 7^ = 7°(cr'^ (g Ts) and for the fourth direction we use 7^. 
We introduce the notations \&t(f)7° = "^(f), H = f d^rH{r) = Y.kH{K) : 



H = Ef=i^^^5o(^(r)7^'^(f + aj) - ^(f + aj)7^'^(f)) + A5o(^(r)7^^(f + 04) - ^(f + 
a4)7^^(^)) + 

M*'(r)*(f) - tEf=iC*(^+ ai)^(r) +^(r)^(r + a^)) 

D D 

H{K) =^{K)[\so^l^sin{Kja) + {M - 2t^cos{Kja))]m{K) (5) 
The eigenvalue of this Hamiltonian are given by: 

e{K) = ±^Yl^^^{2XsoSin{Kja)y + [M -2t ^j^i cos{Kja)Y where D is the spatial 
dimension. For K ^ the system has a gap given by M — 2tD which vanishes at the critical 
value gfj) = 1 and for K ^ [vr, tt, tt, tt] the gap occurs at M + 2tD . In order to investigate 
the effect of the interface we will assume that in the D direction (which we take discrete to 
be j) we have a defect region with the parameters M(j) ,t(j) and \so{j)- For the rest of 
D — 1 dimensions we have a perfect lattice which obeys periodic boundary conditions. We 
note that the eigenfunction for such a problem takes the form Q^{f,j) = e^^'^Us{j) where 
Us{j) is a four component spinor, K and r are the momentum and coordinates in — 1 = 3 
dimensions. The spinor Us{j) obeys the zero eigenvalue equation: 



3 3 
[J2^'stn{K,a) + M{j) - 2t{j) cos(ir,a)]f/,(j) + [t{jh' + XsoimsU + 1) - 

1=1 i=l 

[t{j - 1)7' - XsoU - m^ij - 1) = (6) 

We consider a situation where M(j) = 2t{j) Yl^=i cos{Kia) + e(j) where e(j) = 5j,-ieo — 
5j_ieo. To simphfy the problem we take and \so{3) to be constant. The spinor Us{i) 
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obeys 7^?7±(i) = ±U±{j). Since the zero mode must be normalized we expect that only 
one of the spinors U±{j) will satisfy the zero mode equation . Therefore we have a chiral 
massless fermion at the interface. This condition can be used for any time reversal points in 
the Brillouin zone. For D = 3 the 8 time reversal points are given Ti=ni,n2,n3 = '^in^i, ^2, n-^) 
where rii = 0,1 ,i = 1,2,3 and the gaps occur at K = 7r[?2i, ^2, Tia]. Using the the critical 
condition M = 2tY^^=iCos{Kia) we find that for a fixed value of || the chiral modes will 
occur at K = 7i[ni,n2,n^] only if || obeys || < [cos{it ■ rii) + cos{7i ■ n2)cos{7T ■ n-^)]. Following 
0] we identify the value Si(^ni,n2,n3) = ~1 with the points where ^ < [cos{7r -ni) + €03(71- ■n2) + 
cos{7i-n3)] and Si(^ni,n2,n3) = 1 with the points where ^ > [cos{n-ni) + cos{TT-n2) + cos{7i-n3)] 
. This allows us to make contact with the condition for a strong topological insulator 
:(— 1)'^" = Y[^=i^i{ni,n2,n3), whcu the product is —1 we have a strong TI. Our condition is 
applicable to any interface between two domains, all we need to consider is the equation 
M(j) — 2t{j)[cos{Tr ■ rii) + cos^ii ■ ^2) + cos^tc ■ n^)]. 

The application of these results to D = 3 will allow to treat the D — 1 = 2 boundary as a 



two dimensional surface which gives rise to a Helical metal 32|]. (The particles {e{K) > 0) 
have an opposite chirality in comparison with the antiparticles (the holes). For a TI in 
D = 2 dimensions the boundary corresponds to one dimensional chiral fermions where spin 
up fermions propagate in one direction and spin down electron propagate in the opposite 



direction 



3l|. 



IIB-The formulation of the TI in terms of the Chern numbers 



In order to compute the response of a TI to an external Electromagnetic EM field we 
need to know the non-zero Chern numbers which are the coeficients of the the Chern-Simons 
action obtained after integrating the Fermions. For the remaining part we will consider the 
TI situation with the inverted band condition M—2tD < at i^' ^ . Using the eigen-spinor 
we will obtain the transformed coordinate representation, since this can not be done with a 
single transformation (due to singular points on the torus) we find that the coordinates do 
not commute and therefore the system has a non zero curvature. 

The spectrum of the TI has electron like spinors u{K,s) s=l,2 (positive energy) 
with eigenvalue e{K) > and the antiparticles (hole like) spinors v{K,s) s = 1,2 
for e{K) < 0. This allows to introduce the electron like operator c{K, s) and the 
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— * 

anti — particle (hole) operator b{K, s) which annihilate the fully occupied ground state 
\Go > , c{K,s)\Go >= 0,b{K , s)\Go >= 0. The four component spinor operator \l/(r) = 
[i/{r)i^-^, ^(r)i^|, \E'(r)i^;, \[^(r)i^^]-^ has the momentum expansion: 



^(f) = /" V [c(K, s)uiK, s)e'^ '' + h\-K, s)v{-K, s)e' 



(7) 



The ground state eigenspinor is given by \v{K^s) >. Following we find that 

the transformation \K,a =t, i>— ?■ \v{K,s) > replaces the coordinate representation r* = 
idxi — )■ -R* = idxi + where A~(K) is the connection matrix 

A'{K)=J2 J2 <'v{K,s)\dv{K,s') > \v{K,s) ><v{K,s)\ = 

s=l,2s'=l,2 

D 

<v{K,s)\dKMK,s') > \v{K,s) ><v{K,s')\dK' = 

i=l s=l,2s'=l,2 

Yl l^''A-'^''\K)]dK^ (8) 

i=l o=l,2,3 

The connection matrix A^{K)i forms an SU(2) group. The transformed coordinates 
obey new commutation rules given by the curvature Fij: 

[R\K),R^{K)] = {~i)F,, (9) 
where Fi^j = Ea=i,2,3 and 

i^-^. = 9^.A7(")(i?) - d^,Ai^^\K) + e^,p,A-^^\K)A-^^\K) (10) 



Following 261] we define the n'th Chern number using the notation: Trace[F] = Tr[F] for 



the curvature Fij{K). 

Using the Heisenberg equation of motion for the D dimensional TI we find that the Dirac 



matrices 7 are proportional to the velocity operator '^^jf^^ =7* S2|. In the presence of 
weak external fields an{f,t) we find with the help of the non-commuting coordinates and 
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the relation '^^Jf^'' = Y that the response to the external fields (in agreement with 12] is 
given by: 



D 

c.s. 



where 



(12) 



C.S. 



dt€^ + (9^2 '3'^3 • 



a) The second Ch2 number for TI 



29 



(13) 



43 



44| to derive the second 



We propose to use the theory of first class constraints 
Chern number for the TI. For D = 4 the only non zero Chern number is C/12 
index theorem 26|] relates the nummber of chiral fermions to the Chern number : 



26j. The 



N,-N_ 



(14) 



which allows us to make the connection between the number of zero modes of opposite 
chirality to the Ch^^o. When M—2tD < at -ft" = one finds on the boundary only a single 
zero mode. According to the index theorem this will imply that the Chern number must be 
one. The time reversal points give rise to two connections matrices A~{K) and A~{—K), 
which are defined on the upper and lower side of the torus respectively. Each half torus 
can be mapped to a sphere S^. Using homotopy we have the map 5^ — )■ SU{2). According 



to 



26| we obtain n4[5'[/(2)] = Z2. We conclude that the Chern number Chn=2 = ^2- To 



compute the Chern number we will use A (K) for the North s phe re (upper torus) and 
A~{—K) for the South sphere (lower part of the torus). Following 26| we have the identity: 



Tr[(F)^] = d{Tr[AdA + ^A^]) (F)^ is closed but not exact (we can not satisfy the identity 
with one connection A on the entire torus ). Therefore the C/12 number will be given by a 
difference of two three form Chern-Simons terms: 



Tr[{Ff] = dTr[AdA+ -A^] 

3 

Using the connections A~{K) and A^{—K) we find: 



(15) 
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Jtd=4 21 2tT Js^-Nor-th 2\ iTX Js^-South 2\ iVT 

^ ^ {Tr[{A-^^'''-'^{K)f]-Tr[{A~'^°'''^\k)f]) 



^ ljTr[{A-{K)f]-Tr[{A-{-K)f] = l^, j ^iig-' dgf] = ^2 (16) 

The difference between the two connections A~{K) — A''{—K) = g^^dg defines the tran- 
sition function. At the TRI points [3] the time reversal operator T = —i(T2 ® IK defines 
the matrix elements < ~K, s\T\K, s' >= gs^s'{K) = —gs',s{~K). This relation imposes a 
constraint on the connection field. The solution of this constraint follows from the BRST 



Cohomology 29j. The solution of the constraint will be that g ^dg is an element of the Z2 



group therefore 2^ /^a [(^'^^'^5')'^] takes two values or 1. 
b)The first Chern numbers for D=2 

Following the theory of the Chern numbers for D = 2 the first Chern number is non- 
zero if time reversal symmetry is violated. Therefore the proposal jsO] that the Spin Hall 
conductance is quantized for the two dimensional heterostructure CdTe/ HgTe/CdTe needs 
clarifications. 

bl)- Chi in the absence of time reversal symmetry for D=2 

We consider the Hamiltonian: 

h{K) = KV2 - irVs + M{K)a^ (17) 

When M{—K) = M{K) the time reversal symmetry is broken ( the TR operator is given 
by T = ® o"^)K , where K is the conjugation operator). 

For this case we have a single eigenstate for the ground state \v{K) >. Using this 



eigenstate we compute the connection one form 26| A (K) =< v{K)\dv{K) >. Due to 



the singularity at = 26| we map the upper torus to the North sphere 5*^ defining the 
north connection A~^^°^^'^\K) and the lower part of the torus is mapped to the south 
sphere with the South connection A~^^°^^'^\K). The North eigenvector is related to the 
South eigenvector through a gauge transformation. When we rotate the eigenstate \v{K) > 
by 271 the state is modified by e*'^ ( note the factor ^ which is due to the spin half system). 
Therefore we have the relation between the north and south state : 
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\v{K) >North= \v{K) >South (18) 

which gives: 



— i - = — (f-d^ = - (19) 

2 2 

This imphes that the Hall conductance for the Dirac equation is |^ instead of ^. In the 
presence of a boundary which restricts the system to a; > the term M{K)(j^ will give rise 
to a single spin Polarized massless fermion. 

b2)-The first Chern number for time reversal system when D=2 



Following [30| and the experimental discovery [18|, ll9| the Hamiltonian for the two di- 
mensional heterostructure CdTe/ HgTe/CdTe with the thickness dq^w. exciding 6.3 nm has 
inverted bands and is given by : 



h{K) = KW^ (g) / - K^T^ (g) (T^ + [M - 2t{cos{Kia) - 2tcos{K2a)]T^ ® I (20) 

For ^ ^ the condition for TI is |j < 1. This model has the special feature that for 
spin up and spin down we have two Hamiltonians: 



= (K) = K^T^ - K^T^ + [M - 2t{cos{Kia) - 2tcos{K2a)]T^ (21) 

= (K) = K^T^ + K^T^ + [M - 2t{cos{Kia) - 2tcos{K2a)]T^ (22) 

Therefore both Hamiltonians can be diagonalized separately in the band space r. In each 
band the time reversal operator (does not include the spin operator ) is given by T,- = K and 
obeys = 1. Since and h-^ (separately ) are not time reversal invariant. We conclude 
that in the r space each Hamiltonian has Chi number denoted by Ch[ and Ch\. Following 
the discussion presented above above we have: 
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[RiK),R^iK)] = fl^T'®a' (23) 
Therefore the Chern number will be an integer for each spin polarization: 

(24) 

Contrary to the previous case (61) where the transition function between the North and 
the South pole was \d(f), for the present case the transition function will be dcj) (not half)! 
Repeating the calculation for Ch\ we find Ch\ = —1. Therefore the spin Hall conductance 
will be given by 2 ■ e^/i characterized by the first spin Chern number C/if*" = l[Chl - Chi]. 

IIC-Time Reversal Invariant Topological Insulators in 3 + 1 dimensions in- 
duced by electron-electron interactions reproduce effectively the theory in 4 + 1 
dimensions 

The external fields are restricted to three dimensions. Using the effective action 5Seff. = 
Ch2T^^^ we will take the fourth dimension field 0^=4(7^, t) arbitrarily. We compactify the 
fourth dimension 0(x = X"*^, x^, x^, t) = ^ dx'^a^j_=i{r = x\ x^, x^, x'^, t) [3| and obtain : 

5Sf)j =CK=2- (t){x,t) j d^x j dte^'''''''^^df,a^{x,t)d„ar{x,t) (25) 

There are only two values 6 = Chn=2 ■ <P{S,t) which respect the time reversal symmetry, 
6 = 7T corresponds to a TI and 6 = correspond to a regular insulators. 

We introduce a construction which is entirely in three dimensions space and 
is generated by the electron-electron interactions 

The Hubbard- Stratonovich transformation allows the representation : 

5^p,(f)^(f)7'^v]>(f) + (p,(f))-— i^(p,(fO) (26) 

2K[r,r) 

When the channel P5{r)'-f^ dominates the two body interaction we can use single variable 
q = p^{r). We obtain a four dimensional model build from the three dimensional momentum 
Ki, K2, with the new field (due to interactions) q. This allows us to use the D = A 
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formulation with the vector K = [Ki, K2, K^, q]. (For an arbitrary interaction we can 
construct a four dimensional vector in the following way: 



K = [K,+p,,K2 + p2,K3 + P3,p,] (27) 
Following the procedure from the previous section we compute the second Chern number: 

1 r rgmax. 

= ^ / d^K / dge*'^'*=''Trace[F,,,F,,,] (28) 
Due to the integration measure {paif)) ^ (■^p) (Pa(^0) approximated the range q = p^if) 

to qmin. q qmax.- 

In order to compute the polarization we must consider the external electromagnetic in- 
teractions h^^'^: 

^EM ^ ^^^^^ t) + I ■ M^'"^ (r, t) + (f, t) ■ S (29) 

For a weak electric and magnetic field in the the z direction we introduce E^^^ {Q) = 



(Q) _J Aext.(n\ _ riBI-'-.iQ2,..) 



and A^^^^iQ) 



iQz ^ ' iQ2 

We substitute the transformed coordinates r* — )■ i?* The derivative of the coordinate is 

related to the commutator ^ = El^^\R\R^] ,i = 1,2. We substitute the expression 

for ^ into the magnetic interaction |[^ ■ 6A^^\f,t) + 6A^^'^\r,t) ■ ^] using a slowly 
varying magnetic field, A2^^'{Q) = ''^^^^g^'^^' • Keeping terms to first order in Q we find 
that the interaction can be written as the nonlinear polarizability similar to the axion 
electrodynamics [ll, 12, 60]. 



32 



— dq I d^Ke^'^'^'^TriK [IT' , R^][R\ R'] >)E^{Q)B^{Q) = / d^x / dtE-B 

(30) 



where the expectation value < [R^ , R^][R'' , R''] > is taken with respect to the interaction 
measure Y.aiPair))^^^^iPa{r')) ■ 

III- TRANSPORT ON THE BOUNDARY OF A TOPOLOGICAL INSULA- 
TOR 
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The three dimensional Topological Insulators Bi2Se3,Bi2Tes ,Sb2Te3 support on their 
two dimensional boundary (surface) chiral helical metals described by a two component 
spinor. When the space vector r is rotated by an angle the transformation f ^ r' 
generates a spinor transformation ^(r) — t- \E''(f') = e*'^'^^(r). This implies that when we 
perform a rotation of = 27r the spinor accumulates a Berry phase of vr [9]. This effect 
suggests that in the presence of disorder the state \E'(r) is scattered to a new state \E''(r^). 
When the scattered state is returned back to the point r it will have an extra Berry phase 
of vr. Since localization effects are determined by the interference between the closed path 
and a time reversal path, the extra phase of vr will give rise to destructive interference, 
which means absence of localization. This idea has not been confirmed experimentally. 
The main difficulty being the fact that the experiments use nontrivial geometries such as 
cylinders or spheres. The topology affects the Dirac equation by adding new unknown Berry 
)hases. The Aharonov -Bohm oscillations experiment performed in nano- wires Bi2Se^ in 
35[ seems to support the idea that the maximum conductance should occur for a zero or 
integer flux $ = (in zero flux we have a Berry phase of vr) and the minimum should be at 
$ = (n + where localization is supposed to take place. The experiment of 35| has been 
performed in a cylindrical geometry. On a cylinder the coordinate transformation (from the 
plan to the cylinder) modifies the boundary conditions for the spinors from periodic to 



anti-periodic boundary conditions 



36| . As a result the Berry phase is modified by an 



extra half unit flux. Therefore, when an external flux pierces the cylinder we obtain gapless 
excitation which propagate along the z only when the flux is exactly This means that 
the maximum of the conductance will occur for fluxes $ = (^ + |)^ and the minimum 
should be for zero flux or $ = {n)^. 

Applying an external magnetic field to the two dimensional surface reveals that the 
Landau quantization and the Hall conductance are modified on the boundary of TI. Due to 
the Dirac form the Hall conductance is given by cr^^y^"''^^ = "I" |)- Recent experiments 
of scanning tunneling spectroscopy (STS) [35| reveal the existence of a state at n = which 
is not affected by the value of the magnetic field (the Landau levels for n > 1 scale with the 
magnetic field \/B). The experiments are performed on curved geometries 35] and therefore 
the quantization rules are modified. For example the Hall conductance for the Dirac equation 
on a sphere is a'^'P'^^'^'^{x,y) = 2 ■ a^'^^f"-'^'^{x,y) = (2n + 1)^ [36]. This suggests that for a 
thin sample ct*'**"'(x, y) = (j*°^(x, ?/) + (T^°**°™(x, y) (the measurement might show contribution 
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from the top and from the bottom surface). This suggests that for a thin sample the Hall 
conductance will be similar to the Hall conductance on a sphere. 

The momentum-resolved Landau spectroscopy of the Dirac surface state in Bi2Se3 35 1 
reveals in addition to the Landau spectroscopy, triangular-shaped defects which are caused 
by the presence of vacancies of Bi defects at the Se sites. Clearly the scanning tunneling 
microscopy indicates a modification of the local densities. 

These experiments suggest that the Dirac equation is sensitive to defects and geometry 
and therefore modified quantization rules are expected and probably revealed through the 
new Berry indices 7 7^ |. (In a recent Shubnikov-de Haas oscillations experiment performed 
on the topological Insulator Bi2Se^ the authors present data which show magnetic oscilla- 



tions which correspond to the Landau level quantization En = VF\/'2{n + 'y)heB where 7 is 
the Berry phase which is different from 7 = i.) The appearance of the new indices 77^^ 
might be due to topology of the defects or /and interactions. 

The transport results for the conductance of the TI Bii^j-Shj. show interesting oscillations 
in the presence of a magnetic field as a function of the angle between the axes C3 and C2 
33( 1 . The transport experiments Bii^^Sbx performed in the range 0.7 < x < 0.22 show a 
significant contribution to the conductivity which comes from the bulk. 

It seems that the results are not ruling out the surface conduction contribution, but 
suggests contribution from the bulk. The theory of TI predicts that when a magnetic field 
perpendicular to the two dimensional surface is applied the conductance ax,x decreases 
since the Berry phase which in the absence of magnetic field is vr, becomes in a magnetic 
field 7r(l ± , ^ =). Therefore backscattering is allowed ( The parameter A contains 
two contributions, the normal component of the external magnetic field and the magnetic 
exchange which is due to the magnetic impurities. Both effects can be described by a 
modification of the Hamiltonian by 6H = A ■ a^. It is known that in a system with spin 
orbit interaction the conductivity also decreases in the presence of a magnetic field j42|. 
Therefore a quantitative analysis is needed in order to separate the two effects. If the 
Topological Insulators are to be responsible for the conductivity one needs to understand 
the surface-bulk contributions. 



IIIA- The Effect of Static Disorder on the Second Chern Number 



16 



The transport experiments 33|] performed on the mixed crystal Bii^^Shx TI for 0.7 < 



X < 0.22 show that a significant contribution to the conductivity comes from the bulk. In the 
absence of disorder the condition < 1 guarantees the existence of the TI. The presence 
of disorder, is included trough the impurity potential: Vdisorder = Vo(r)7° + Vi{r)I (Vo^r) 
represents the space dependent mass term - the orbital splitting energy) . The mass gap M 
is normalized by both terms. 

a) -The two mode crystal 

For inhomogeneus disorder we have a two mode crystal, namely the crystals splits into 
domains where the mass gap M takes values, ^^^^^ < 1 (domain^) and > 1 

(domain C). At the interface between the domains we will have chiral massless fermion 



propagation. The recent transport experiment |33| hints that a significant contribution to 
the conductivity comes from the bulk. Therefore the two mode proposal might explain the 
high conductivity coming from the domain walls in the bulk. 

b)-The homogeneous virtual crystal. 

When the system is homogeneus disordered, the value of the mass gap can be replaced by 



y/ < NP > = M +plus Gaussian fluctuations of M — M. The value ^ = 1 represents the 
transition point between the TI and a regular insulator. Due to the Gaussian fluctuations we 
expect the following transitions : The region < 1— e is a TI. The region 1 — e < < 1+e 
corresponds to a chiral metal which will become a regular insulator for > 1 + e. 

In order to perform a quantitative analysis we will consider the Hamiltonian Hdisorder = 
H + Vdisorder- Usiug the cigen-spiuors \v{K,s) >,\{uK,s) > we compute the transformed 
disorder matrix Vs,s'{K, K') (with the matrix elements vfgf{K,K'), , vfg'!'{—K,—K'), 
vff\K, —K'), vf^f\—K,K') computed with the help of the Fourier transform of the 
disorder potential and the eigen-spinors ). 

H = [ Y.'^e{K){c\K,s)c{K,s) + b\K,s)b{K,s)] (31) 
Vdisorder = f Y^f (cHK , s) , biK , s))[Vs,siK , K')MK' , s]) , b^K' , s')f (32) 
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We will compute the self energy and we will extract the effective mass gap M. Pertur- 
bation theory can be used since time reversal symmetry guarantees that the matrix elements 
Vs^s'iK, K' = —K) are zero. In the presence of disorder the ground state spinor \v{K, s) > 
is replaced by \\{K,s) >: 



s'=l,2 p ^\^) ^\P) 



l\\2 



^ \V{K,P)V*{K,S)U{P,S)\ 

sh^P <K)-e{P) 
\v{K,s) > + ^ ^S(^'^')(^,POI^(^',«') > 

s'=l,2 p 

+ E E^^''''^(^'^')k(^',«')> (33) 

s'=l,2 p 

In order to test the possible phase transition for the Chern number, wc propose to compute 
the new curvature F which replaces F. For this purpose we introduce the coordinate X* = 
idxi and compute the new connection: Aj~. 



A- =< v(i^,si)|X*|(v(i^,S2) >=< v{K,s^)\X'\v{K,S2 > + 

E E Y^{b^'''''\k,p))*x'b^'^^'"\k,p)) 

s'=l,2s"=l,2 p 

+ E E Y^{c^'''''\K,p)yx'c^'^''"\K,p)... 

s'=l,2s"=l,2 p 



(34) 



Using this representation we obtain the curvature matrix: Fj^- = [X-' , ]. 
Using the new curvature Fj j we obtain the second Chern number: 

C/l2 — / Tr < [777(7") ] > average— disorder (35) 

In order to answer the question if the clastic disorder destroys the TI we will perform 
analytical and numerical calculation. In order to modify the Chern number by disorder 
the representation for the new eigenvector |v(i^, s) > must be divergent! When the series 
converges we find that [X^, A~] — [X^, A~{K)i\ and the second Chern number is unchanged. 
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In order to determine the transport properties we will construct an effective non-linear 
sigma model using the Keldish formalism j48 |. 



III.B. The TI with Disorder Using the Keldish Method 

In order to compute the conductivity we need to construct an effective non-linear sigma 



model using the Keldish formalism. Following [48|| we introduce two fermionic fields which 
reside on the forward and backward contour. We replace the four component spinor \l/(r, u) 
by a eight component spinor ,^(r, w) = [^(r, w), \l'^(r', — o;)]"^. Using this formalism we 
represent the disorder action S disorder'- 

S,,.„,„ = / d-rj d.V._,..,(n(*'(F,u,)E3*(f..) (36) 

Performing the average over disorder we can introduce a set of Hubbard Stratonovich 
fields [50|, l58| Qjif, {uj,uj') which act as gauge fields in the space S3 ® 7-' (S3 acts on the 
frequency space). Using the Hamiltonian given in section II A (the TI in D dimensions ) + 
disorder we have: 



S= d^r dco dco'l J2 *^(r,t^)(S3®7^)((-i)4- 
J J i=i,..-D 

+Qj{f, oj, oj')^{f, tu') + a;)(S3 ® -f^)i{M + Qo{f, oj, Lo')^if, to') 

-t *^(r,u;')(S3®7')V,^*(r-',a;')] (37) 

j=l,..D 

where V'j'i/{f,u) = [^(r + aj,ijj) + ^(r — aj,uj) — 2^(r — aj,uj)] and the Gaussian gauge 
fields are controlled by: Hq = f d^r f du f du'U2j=i5 2^{Qj{i^,{^,<^'))'^]- 

This effective action is used to construct the effective action for transport. Integrating the 
fermions we obtain the effective action in terms of the gauge fields Q^. The effective action 
is a function of the curvature F^^^{co, u') = dyQ^if, u, uj') — d^Q^if, u, co') and Chern-Simon 
term. 

Using the Renormalization Group we will establish the phase diagram: In the metallic 
phase the Maxwell term vanishes . If the coefficient of the chern number obeys, C/;,„=2 7^ 
we have a chiralic metal. When Chn=2 = we have a regular metal. When the Maxwell 
term is non zero we have either a regular insulator Chn=2 = or a T.I for Chn=2 7^ 0. 
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III.C. Interaction and Disorder 



For the one dimensional boundary TR symmetry forces the backscattering term to vanish. 
Since umklapp is allowed jsil the presence of the Coulomb interactions for D = 1 drives 
the system to the localized phase. We have found that the Luttinger parameters decreases 
511 . 52I causing the chiral edge excitations to localize. The two dimensional surface 



below T 



represents the boundary of the three dimensional Topological Insulator which is described 
as D = 2 chiral metal given by h{K) = v{K'^ai — K^(y2) 

where K is the momentum and a are the Pauli matrices for the spin of the electron. This 
Hamiltonian has two eigen-spinors \u{K) >= [u^{K) , ui{K)]'^ for particles and \v{K) >= 
[v^{K),vi{K)]'^ for antiparticles which have opposite chirality. The two component spinor 
operator is given by: \E'(r) = [${r)f,'${r)i]'^. Using the ground state \Go > we introduce 
the particles and anti-particles operators: c{K)\Go >= 0, b{K)\Go >= 0. For this case we 
observe that the spinor operator is given by the chiral Weyl form: 
^(f) = J^,[c{K)u{K)e'^-'' + b^{-K)v{~K)e'^-'] . 

The experiments show that the Dirac Cone is below the conduction band. This can be 
considered by introducing a positive chemical potential fi > 0. The presence of the chemical 
potential changes the scaling dimensions of the operators and can cause the impurity and 
interaction potentials to become relevant variables (in the sense of the Renormalization 
Group analysis). At the edge points of the Fermi sphere kp = [0, ±|fci;'|] , kp = [±\kp\,0] 
the electronic spectrum is quadratic! (The excitations in the vicinity of the Fermi surface 
are given by e{q) = ■ q + r^f-r) where g is a measure with respect the Fermi energy.) 



HID- Topological Insulators in Curved Spaces 



It seems that in curved spaces caused by geometry or topological disorder such as dislo- 
cations and disclinations we can have the following scenario. For a TI we dislocations might 
create boundaries such that chiral excitations might be bound to the dislocations. On the 
other hand if we have already massless fermions some of the topological defects might cause 
localization. 



In order to investigate the effect of the 
a geometric formulation for curved spaces 



bpological Disorder on the TI we will introduce 



36 



37 



37|, 



54 



55 



-57 



6l| . We consider an n 



dimensional manifold where a mapping from the curved space X", a = l,2...n to the local 
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flat space exists, x^, = 1, 2...n. We introduce the tangent vector 37| e°(x) = ^^J^^ which 
satisfies the orthonormahty relation e^(x)e||^(x) = Sa^h (here we use the convention that we 
sum over indices which appear twice). Representing the curved coordinate in terms of local 
fiat coordinates gives rise to the metric tensor e^(af)e"(x) = g^^v{x)- We introduce the linear 
connection which is determined by the Christoffel tensor ^ which is defined by: 



= -^iudx (38) 

Next we introduce the vector field V = V^da = V^df^ and the spinor \l'(x) . The covariant 
derivative of the vector field V"' is determined by the spin connection u^f 

D^V\x) = d,V\x) + u'^^.V'' (39) 

Similarly in the basis x^ the covariant derivative D^V^{x) is determined by the Christoffel 
tensor: 

D,V''{x) = d,V''{x) + Tl,V^ (40) 

The spinor connection determines the covariant derivative D^V°'{x) defined above 
and as well determine the covariant derivative of the spinor '^{x) given in term of the gamma 
matrices 7. 

D,^>{x) = {d, + \ufYl'^,^'])^{x) (41) 
We first evaluate the Christoffel tensor which is a function of the metric tensor g^ i,{x^. 

^iu = -\ Yl 9^''i^M9uA^) + 9^9uA^)-dr9^.A^)] (42) 

T=l,2,..n 

The relation between the spin connection and the linear connection can be obtained from 
the fact that the two covariant derivative of the vector V are equivalent. 

D^V' = e^D^V (43) 
Since we have the relation = e^V^ it follows from the last equation 

D^[el] = D.d^e'' = {D,d,)e^ + 9,(D^e'^) = (44) 
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Using the definition of tlie Cliristoffel index we find the relation between the spin con- 
nection and the hnear connection: 



D,[e^J\ = d,el{x) - T^XM + a;^,,et(x) ^ (45) 
From this equation we determine the spin connection: 



_leA-e'^'''(ape.,,-9<,ep,,)e^ (46) 
This allows to rewrite the Dirac equation h = '~f°'{—ida) as : 

^curved ^ . Y^_,^^d^ + ^^a»l^^^ ^,]) (47) 

For the mapping of a surface to another curved surface the transformation is simplified 
. For the surface case n = 2 the square of the length element [dx'^ + dy'^] on a flat surface 
is related to the square length element ds"^ for a curved surface in the following way: ds"^ = 
e~^^^\dx'^ _l_ g^y2j ^]2grg is the conformal transformation from the flat to curve surface. 



36| to compute the form of the Dirac equation for 



This method ha been used in reference 
a sphere. 

c)The effect of Topological disorder on the Topological Insulators- Disloca- 
tions and Disclinations-Computation of the connection. 

The presence of dislocations and disclinations deform substantially the lattice structure 



and therefore might affect the electronic transport. Following |39l] the deformation of the 
crystal give rise to a coordinate transformation which will affect the Dirac equation. Given a 
perfect lattice with coordinates x*, lattice deformation will modify the coordinate to X"(x) = 
x"" + e'^(x) where e"(x) is the local lattice deformation. As a result the metric is given by 

n (7\ - Q^" 

Screw dislocations are characterized by the Burger vector B which is parallel to the 
dislocation line contrary to the edge dislocation where the Burger vector is perpendicular 
to the dislocation line . The metric tensor for a screw dislocation is : ds"^ = (dz + ^dOY + 
dp'^ + p'^dcf)'^. 
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Using the tangent vector we represent the burger vector in term of the torsion: 

jdX^el ^ J J dX^'dX^id^et - d^e^) = -B^"'^ (48) 

In order to describe disdination which has no torsion but has curvature. The metric 
tensor for a dischnation is given by ds^ — dz^ + dp^ + a^p^dcf)^ The effect of dischnations 
can be understood as an angle deficit 27r • a ,therefore once a rotation of 27r is performed the 
spinor will rotate by an angle tt • a < tt . Therefore we obtain a Berry phase of tt • a which 
will allow backscattering and will allow localization. ( Quantitatively the effect of the line 
deficit can be expressed as the strength of a vortex.) We believe that the dislocation might 
play an important role for computing transport properties, contrary to regular disorder some 
of the dislocations and dischnations can cause localization. 

d)An explicit demonstration of the method presented for an edge dislocation 
in two dimensions 

We consider the two dimensional chiral Dirac equation in the absence of any dislocation: 

h^a\-id2) + (T\idi) (49) 

We consider an edge dislocation with the burger vector B^'^'> ( the value of the vector is 
the lattice constant in the y direction). 

The dislocation is given in terms of the torsion tensor T^=i^,^=2 

T;zI=2 = d,el - d,e\ = S^5^(r) (50) 
Using the solution of the Laplace equation we obtain the tangent components: 

^l-^T^^, ^2 = 1-/3^T^ (51) 
[x-^ + y^) [x^ + y^) 

where /3 = and e]; = 1 ,62 = . 

Using the tangent components we obtain the metric tensor to first order in j3: 

gii = 1, 512 = /3 2 ^ 2 ' 922 = 1- /3 ^ , 921 = (52) 

j^yz x'^ -\- 

Using the the vectors components and the results obtained in previous equations, we 
compute the transformed Pauli matrices: 
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a\x)^a\ a\x) ^ - f3^^y (53) 
The spin connection is obtained from equation 46 — 47: 

Col^ = -ool^ ^ TT/SS^X) (54) 

As a result the chiral Dirac equation in the presence of an edge dislocations becomes: 

o X I o 

Using the new Dirac Hamiltonian in the presence of an edge dislocation we obtain the 
following eigenvalue equation for the two component spinor ,E is the energy , /i is the 
chemical potential and f3 = ^(jedge- 

{E - /x)t/t(f) = [(a. - idy) - e-'fg,,gJ\r)]U^{r) - ^^^,,^,(-^)9.f/^(r) , 

The presence of the complex coupling constants e^cjedge restricts the solutions. Preliminary 
results indicate that one dimensional propagating solution along the line y = —x exists, and 
in addition we The interesting one dimensional solutions confined to the dislocation line 
is interesting since it shows that the conductance will be dominated by one dimensional 
channels confined to the edge dislocation. 

Therefore for a number of dislocations centered at points i = 1, ■■■N we will have 

N one- dimensional conducting channels confined to the line {y{x) — yi) = —{x — Xi). 

e)The Chiral metallic state induced in a Topological insulator by a Dislocation 

Using the coordinate transformation induced by the dislocations we can rewrite the 
model given in section II A. The new model will contain the Dirac equation in 4+1 dimen- 
sions curved space (due to dislocations) plus the non relativistic mass term, M'^ {r)^> {r) — 
^X]jLi(^(^+ 'j^i)^(r) + ^(^)^(^+ CLj)) which in the continuum will be replaced by the 
Laplacean = A which in the curved space will become the Laplace-Beltrami operator 
^LB 0, Is^]. The dislocation is believed to give rise to a chiral zero mode (along the 
dislocation line). The eigenvalue equation for the topological insulator in the presence of a 
dislocation can reveal this fact which has to be further investigated. 
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IV. Topological Z2 Realization for the + %'Py Chiral Superconductor Wire 
Coupled to Metallic Rings in an External Flux 

The experimental realization of the "p^ + Wy pairing order parameter (p-wave supercon- 
ductors) in Sr2RuO/^ , ^He — A and the ^ = \ case represents interesting new excitations. 
In these (weak coupling) + %'Py pairing the excitations are half vortices which are zero 
mode energy Majorana fermions. We consider a p-wave superconductor confined to a one- 
dimensional wire. At the edges of the wire a; = and x = L the pairing order parameter 
vanishes and two zero modes appear at the edges (x=0 and x=L). Due to the charge conjuga- 
;ion of the Bogoliubov spectrum these zero modes should be Zero mode Majorana Fermions 



53|. Mapping the problem to the Ising model one can show that the ground state is a Z2 
doubly degenerate ground state. As a result the single particle excitations are non-local. As 
a concrete example we couple the p-wave wire to two rings which are pierced by 
external fluxes. In the first stage we will study a spinless model and in the second stage 
we will consider the effect of spin. 

The model for the spinless case is .H = Hp^w + Hi + H2 + i7(p_vK,i) + -f^(p-H/,2)- 
Hp-w is the p-wave superconductor with the pairing gap A and polarized fermion 
operator C{x) = C^=^{x). Hp_w = E.,.'(C+(a;)^(a;') + h.c.) + XE,C+{x)C{x) - 
'^J2x x'i'yx,x'C^{x)C^{x') + 'yx^x>C{x')C{x)]. The time reversal and parity symmetry are 
broken, 7^ ^/^^.^^^ = —j^ ^i^^^^. Following [53| one introduces the Majorana fermions 



) V2{x) [53| and finds: C{x) = ?7i(x) -|- irj2{x) ; C^{x) = rii{x) — irj2{x) and 



?7i(x) = |(C(x) + C+(a;)) and rj2{x) = ^{C{x) — C+(x)). The Hp_w Hamiltonian has 
the pairing boundary conditions A(x = 0) = A(x = L) = and A(x) = Ap for 
< X < L. This model is equivalent to the Id-lsing model in a transverse field 53| : 
Hp-w = t / dxfj^lo'idx + zA(x)o-2]57 where if = {■r]i{x) , t]2{x)) . 

The two eigenfunctions are charge conjugated with two bound states at x=0 
and x=L. 

rii{x) = rjL^^{x) + .... and rj2{x) = riQ^^{x) + .... 

$°(x) and $-^(x) are the bound states at x=0 and x=L and 770 and rji are the zero modes 
Majorana operators. They obey (r^o)^ = Vo ; (^0)^ = | and (t^l)^ = rji , (?7l)^ = | .This 
suggests that the ground state is Z2 doubly degenerate. 

Using the zero modes Majorana Fermions rjQ and rji we obtain a fermion operator q = 
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|(^L + ^^o) and = — irjo) which obey anti-commutation relations [q, = 1 which 
obeys g^|0 >= |1 > and q\l >= |0 >. The p-Wave wire is described by the ground state 
|0 > when the wire contains only pairs. When one electron is added to the wire the state is 
|1 > . The energy of this state is given by e . The value of e depends on the length of the 
wire, when L — oo we have e — )■ IsS]. 

Hi + Ho represent the Hamiltonian for the two metallic rings each pierced by a flux 0i 



22|. 



H,= &-^ dx ^4(x)i-^^, 



2lT 



1,2, 



[2m 

In the long wave approximation the p — Wave Hamiltonian Hp^w is replaced by Hp^w 



eq^q . 

H{p-w,i) + -f^(P-w^,2) represents the coupling between the p — Wave wire to the two rings 
trough the matrix element g. 

H(p-w,i) + H(^p^w,2) = igiqKi'i + 'ipi + i'l- ^2) + qii'l + ipi-'il>l + V'2)] 

Using the Heisenberg equation of motion we substitute q{t) and and reduce the 
problem to a eigenvalue problem in terms of the rings spinless electron operators ipiix), 
il'2{x)- We will work with the Fourier transform in the frequency z representation. We 
replace the two electron operators by the four component spinor, 

^ix,z) = ['ipi{x,z),'ip2ix,z),'iljl{x,-z),'iljl{x,-z)f. 



( ^i{x,z) \ 

ij2{x,z) 

i^tix, -z) 
\iljt{x,-z) j 



(56) 



ei 

£2 

-62 

-ei 



+ 



-z e e z 
e —z —z — e 
e —z —z — e 
z —e — e —z 



We use the notation: = {—id^ + ^iY- 
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A full solution involves a numerical computation for the eigenvalues and eigenvectors. 
We find that the problem simplifies in two limits e » g and e « g. 

a) When e « g the current in each ring will depend only on the flux applied. For equal 
fluxes we will find that the current and the magnetization are proportional to the applied 
flux. This case corresponds to the Andreev reflection at each interface. 

b) When e » g the current through the current and the magnetization will be propor- 
tional to the {fluxy. This case corresponds to a crossed Andreev reflection 

We believe that measuring the current and the magnetization as a function of the flux 
will allow an indirect conflrmation for Major ana fermions. 

V-CONCLUSION 

This article aims to introduce a variety of powerful geometrical methods for condensed 
matter systems such as connection and curvature. Electrons in sohds are characterized by 
Bloch bands. When the band crosses one obtains non-trivial phases which gives rise to non 
-commutativity in momentum space characterized by the Chern numbers. The Topological 
insulators are material with spin-orbit interaction. At special point in the Brillouin zone the 
model is time reversal invariant . As a result a new type of insulators characterized by the 
second Chern number is identifled. At the boundary of the Topological Insulators an odd 
number of Chiral states are identifled. Those states are robust againt disorder and represent 
perfect metals. 

The idea of protected boundary can be generalized to any interface including disorder or 
Topological defects such as dislocations and disclinations. 

We introduce Majorana fermions and show that they can be detected by measuring the 
persistent current in metallic rings attached to a P-wave superconductor. 
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